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The spatial distribution of a species provides valuable in¬ 
sights into how it uses its habitat. In general, species are 
found at sites that provide resources adequate for their 
survival and reproduction (Grinnel 1917, Hutchinson 
1957, VAN Buskirk 2005, Bowles et al. 2006, Indemaur 
et al. 2010, Menin et al. 2011). The way in which a spe¬ 
cies uses these resources, such as specific microhabitats, 
may vary considerably in both time and space, resulting 
in a heterogeneous distribution of individuals within the 
environment (Pianka 1973). This variation in the distri¬ 
bution of individuals has been studied widely in all ver¬ 
tebrate groups, including amphibians. In amphibians, 
researchers have observed that environmental factors 
(e.g., the structure of the vegetation and microclimate) 
determine the spatial distribution of individuals (Spot- 
ila 1972, Crawford & Semilitsch 2008, Correa et al. 
2012). These factors may be affected by anthropogenic im¬ 
pacts, such as forest fragmentation leading to edge effects. 
The response of the species to these impacts in terms of 
the distribution of species may be either positive or nega¬ 
tive, or neutral (Murcia 1995, Ries et al. 2004, Harper 
et al. 2005). 

Some species of salamander (Caudata) respond nega¬ 
tively to anthropogenic disturbance (Messere & Ducey 
1998, Harper & Guynn 1999, Marsh & Beckman 2004). 
Salamanders are amphibians that occupy an enormous 
variety of habitats (e.g., bodies of water, forests, caves) and 
microhabitats (e.g., leaf litter, rocks, vegetation) in various 
parts of the world (Vitt & Caldwell 2014). In arbore¬ 
al salamanders, the adults and juveniles typically occupy 
distinct strata in the vegetation (Waldron & Humphries 
2005), with seasonal variations in the use of substrates by 
the different individuals (Sequeira et al. 2001). 


A considerable diversity of salamanders is found in the 
Neotropics, and the genus Bolitoglossa is the most diverse, 
with more than 100 known species (Frost 2015), which are 
found in forest habitats, primarily in the leaf litter and veg¬ 
etation (Duellman 1978, Chan 2003, Ortega et al. 2009, 
Neckel-Oliveira et al. 2011). Despite the overall diversity 
of Neotropical salamanders, most ecological studies have 
focused on reproductive traits or microhabitat use (e.g., 
Houck 1977, Bruce 1997, Chan 2003, Ortega et al. 2009, 
Neckel-Oliveira et al. 2011) and little is known of their 
distribution in the habitat. 

Five species of Bolitoglossa are known to occur in the 
Amazon basin of Brazil, including Bolitoglossa paraensis 
(Brcko et al. 2013). This species is found in the understo¬ 
rey vegetation (Neckel-Oliveira et al. 2011), where its 
ethology is influenced by climatic variables, such as rain¬ 
fall and relative humidity (Correa et al. 2012) . Populations 
of B. paraensis are known to occur in forest fragments in 
the northeast of the Brazilian state of Para (Brcko et al. 
2013), which is an area of intense deforestation due to 
farming and urban growth (Gerwing 2002, Souza et al. 
2003, Borner et al. 2007, Pacheco 2009). We investigated 
the variation in the distribution of Bolitoglossa paraensis 
specimens during different periods of the year. The goal 
of this paper is to answer two questions: (i) is the distri¬ 
bution of specimens affected by the distance to the forest 
edge? And (ii) is the distribution of the specimens related 
to their size? 

We conducted our study in the Gunma Ecological Park 
(CEP), a 540-ha forest fragment in the municipality of San¬ 
ta Barbara do Para (oi°i3’86” S, 48°i7’4i” W) in the north¬ 
east of the state of Para, Brazil. The CEP is covered mainly 
by dense ‘terra-firme’ rainforest (61% of the area), with can- 
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opy trees of more than 25 m in height. There are also areas 
of floodplain (5%) and secondary forest (34%) (Amaral 
et al. 2009). The GEP is intersected by two dirt roads and 
the PA-391 state highway. Perpendicular to the dirt roads, 
there are four trails (< 1 m wide), more than 400 m from 
one another. 

The local climate is humid equatorial, with a marked 
rainy season between December and May, and a period of 
reduced precipitation between June and November (Albu¬ 
querque et al. 2010). During our study period, the Na¬ 
tional Meteorological Institute (INMET) recorded month¬ 
ly precipitation at the nearby town of Belem of between 
55.7 mm in November and 662.2 mm in March (Fig. 1). At 
the GEP, the median temperature inside the forest varied 
from 24.1 to 25.i°C, while the median humidity ranged at 
86-95% during our study period (Gqrrea et al. 2012). 



Months ( 2009 ) 


Figure 1. Monthly rainfall at Belem as obtained from a database 
made available by the Meteorological National Institute (INMET, 
http://www.inmet.gov.br). 


We demarcated twenty-three 250-m^ (5 x 50 m) plots 
within the area of the park, at a distance of at least 200 m 
from each other. This method was adapted from Rodel 
& Ernst (2004), and appears to be the most effective for 
long-term studies of terrestrial amphibians. The plots 
were located randomly within the existing GEP trail sys¬ 
tem and georeferenced with a Garmin 76 WAAS GPS re¬ 
ceiver. We used the georeferenced position of the plots to 
locate them on a schematic map of the study area and cal¬ 
culated the distance of each plot to the nearest forest edge 
(Fig. 2). 

We collected data during four periods in 2009: the mid 
rainy season (March-April), late rainy season (May-June), 
mid dry season (August), and late dry season (October). 
On some nights during the rainy season, the rainfall was so 
heavy that it impeded fieldwork. For this reason, two plots 
were sampled during the first two days of April (mid rainy 
season) and three days during the first three days of June 
(late rainy season). While June coincides with the onset of 
the dry season, we considered the data from this month 
(taken on the first three days) to be representative of the 
late rainy season. 

We surveyed all 23 plots once during each period for 
the presence of Bolitoglossa paraensis, based on active vis¬ 
ual searches conducted between 19:00 and 22:00 h. During 
each survey, three people systematically walked through 
the plot looking for salamanders on the different sub¬ 
strates known to be used by this species, such as the sur¬ 
faces of leaves, herbaceous plants, lianas and tree trunks up 
to a height of 3 m (Neckel-Oliveira et al. 2011). When 
an individual was found, we recorded the period of the 
year, plot, height above the ground, and snout-vent length 
(SVL). We measured the height above the ground for each 



Figure 2. Map of the Gunma Ecological Park in northeastern Para, Brazil, showing the distribution of the sample plots. Legend: grey - 
forested; white - open areas; black marks - sample plots; dashed lines - highways and roads. 
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individual with a 5-m measuring tape and used a 0.1-mm 
digital calliper to measure SVL. 

We analysed our data separately for each of the four 
study periods. Because they did not assume a normal dis¬ 
tribution according to the results of the Shapiro-Wrlk test 
even after transforming them, we analysed our data using 
non-parametric tests, as suggested by Zar (loro). 

As individuals may be distributed homogeneously or 
heterogeneously within the forest, we compared their ob¬ 
served distribution (number of individuals in each plot) 
with the null-expected distribution using a Chi-square 
test. The alternative hypothesis indicates a heterogene¬ 
ous distribution of individuals within the forest. During 
each period, we recorded ‘outliers’, that is, data collected 
in atypical nights (when it rained for up to 30 min prior to 
the surveys), as proposed by Correa et al. {1012). While 
we collected data following the standard procedures dur¬ 
ing these nights, we ran Chi-square tests with and with¬ 
out ‘outliers’ in order to test their reliability. In addition, we 
used a Spearman coefficient to test for correlation between 
the number of salamanders in each plot and their distance 
to the nearest forest edge. 

We tested the variation in the median height above the 
ground at which salamanders were found between the dif¬ 
ferent sample periods using a Kruskal-Wallis nonparamet- 
ric analysis of variance. When the null hypothesis was re¬ 
jected (p < 0.05), we used a multiple comparison test to 
identify which periods contributed signiftcantly to the 
observed variation. We also used a Spearman coefficient 
to test for correlation between the size of individuals and 
their height above the ground. We ran all analyses in PAST 
(Hammer et al. 20or) and the R packages Pgirmess and 
Stats (R Development Core Team 201^). 

In the 23 plots, we recorded 90 individuals of Bolito- 
glossa paraensis in March-April (mid rainy season), 46 in 
May-June (late rainy season), 45 in August (mid dry sea¬ 
son), and 25 in October (late dry season). When the out¬ 
liers were included in the analyses, the individuals were 
distributed heterogeneously over the plots during all four 
periods of the year (Table r). When the outliers were ex¬ 
cluded, however, a heterogeneous distribution was found 
only for the mid rainy season (Table r). No correlation was 
found in any period between the number of salamanders 
per plot and the distance to the forest edge (Table 2). 

The median height above the ground of B. paraensis 
specimens varied signiftcantly over the study period (H = 
39.924; p < o.or), with differences in the mid rainy season 
and the two periods of dry season (p < o.or). We also found 
a negative correlation between salamander size and height 
above the ground in the mid rainy season, whereas no cor¬ 
relation was observed in any other period (Table 2, Fig. 3). 

In tropical amphibians, variation in the distribution of 
specimens is usually associated with the breeding season 
and the availability of prey, which vary according to climat¬ 
ic conditions (Watling & Donnelly 2002, Chan 2003, 
Whitfield & Donnelly, 2006, Rodrigues & Santos- 
CosTA 2014). In this context, some ecological characteris¬ 
tics of Bolitoglossa paraensis, such as reproduction (Neck- 


Table 1. Results of the Chi-square tests for the differential distri¬ 
bution of Bolitoglossa paraensis in the plots in the four periods of 
the year surveyed in the Gunma Ecological Park. 


With outliers Without outliers 

P-value P-value 


March-April 

90.933 

< 0.01 

47.536 

< 0.01 

May-June 

34 

0.049 

29.218 

0.257 

August 

41.378 

< 0.01 

18.657 

0.607 

October 

45.84 

< 0.01 

25 

0.247 


Table 2. Spearman correlation coefficients for the distance to the 
forest edge vs. the number of individuals per plot (DEE) and sala¬ 
mander size vs. height above the ground (BSH) during each of 
the four study periods. 




DEE 

P-value 

r 

s 

BSH 

P-value 

March-April 

-0.118 

0.598 

-0.224 

0.028 

May-June 

-0.316 

0.139 

0.114 

0.436 

August 

0.069 

0.781 

0.011 

0.947 

October 

0.175 

0.419 

0.259 

0.201 


el-Oliveira et al. 2on) and abundance patterns (Correa 
et al. 2012), are influenced by climate in some periods of 
the year, but not in others. In our study, we observed a sim¬ 
ilar pattern in the distribution of salamanders within the 
forest. 

We observed that the distribution of B. paraensis within 
the forest is not homogeneous, at least during mid rainy 
season, when rainfall levels are highest (March-April, 
Fig. r). Plethodontid salamanders tend to use microhab¬ 
itats with high relative humidity to avoid dehydration 
(Spotila 1972). Bolitoglossa paraensis appears to be espe¬ 
cially sensitive to the variation in relative humidity, given 
that more active individuals were observed during periods 
with increased humidity (Correa et al. 2012). 



Figure 3. Correlation between Bolitoglossa paraensis size and the 
height above the ground of the specimens encountered during 
surveys. 
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Bolitoglossa paraensis reproduces during the rainy and 
early dry seasons, when the relative humidity inside the 
forest is relatively high (Neckel-Oliveira et al. 2011). 
During the breeding season, salamanders aggregate at spe¬ 
cific points in the forest, influenced primarily by two fac¬ 
tors: microhabitat quality and the presence of conspecifics, 
as observed by Quinn & Graves (1999) in a population 
of Plethodon cinereus (Plethodontidae). This is typical of 
aggregations associated with the distribution of resources 
in the habitat (Pounds & Crump 1987, Hews 1993, Ber- 
TONA & Chiaraviglio 2003), which may be food, mates 
or breeding sites. 

While the distribution of salamanders within the for¬ 
est was not homogeneous, we did not find any evidence of 
an edge effect. In this case, B. paraensis may be considered 
to be a non-edge-avoiding species in the classification pro¬ 
posed by Lethinen et al. (2003). In the tropics, the effects 
of forest fragmentation have been widely studied in anuran 
amphibians (e.g., Vallan 2000, Lehtinen et al. 2003, 
Dixo & Martins 2008, Correa & Rodrigues 2015), but 
not salamanders, which have only been studied in a tem¬ 
perate region (Gibbs 1998a, b), where some species pre¬ 
sented a negative response, while others were unaffected. 
However, while B. paraensis is unaffected by the distance to 
the edge, it does need forested habitats, in which it is typi¬ 
cally found on the surfaces of the leaves of trees and her¬ 
baceous plants (Neckel-Oliveira et al. 2011). As in oth¬ 
er plethodontid salamanders (Spotila 1972), B. paraensis 
is also affected by relative humidity (Correa et al. 2012), 
which is much higher in forest habitats than in open ar¬ 
eas (Chen et al. 1999). Our results thus differ from the ex¬ 
pected pattern for salamanders, given that temperatures 
are usually higher near the forest edge, while relative hu¬ 
midity will be lower here than inside the forest (Saunders 
et al. 1991, Murcia 1995), and that such changes in micro¬ 
climate are known to have a negative effect on plethodon- 
tids (Spotila 1972), including B. paraensis (Correa et al. 
2012). However, at the CEP, relative humidity does not ap¬ 
pear to vary in relation to the distance from the forest edge, 
at least not enough to affect the distribution of B. paraensis 
within the forest. 

During the mid rainy season (March-April, Fig. 1), the 
larger salamanders tended to be found lower down in the 
forest, in contrast with the positive relationship between 
size and height above the ground observed in Bolitoglossa 
cf pandi (Ri'o-Garci'a et al. 2014), and the lack of any re¬ 
lationship in other Bolitoglossa species (Mead & Boback 
2002, Ortega et al. 2009). In these studies, however, data 
were collected throughout the year, rather than the season¬ 
al analyses conducted here, and these differences may con¬ 
tribute to the differences found in the results. It is impor¬ 
tant to note that the ecological characteristics of amphib¬ 
ians may vary considerably between seasons (Watling & 
Donnelly 2002, Chan 2003, Whitfield & Donnelly 
2006, Correa et al. 2012, Rodrigues & Santos-Costa, 
2014) and that this variation may not be detected if the data 
are not analysed separately. 


The fact that larger salamanders were found more fre¬ 
quently at lower heights during the mid rainy season may 
be accounted for by the B. paraensis breeding season. In 
B. nicefori, gravid females are much heavier and may be 
less able to climb vegetation, and thus remain low down 
in the substrate or on the leaf litter, where they lay their 
eggs (Ortega et al. 2009). This is consistent with the fact 
that B. paraensis reproduces mainly during the mid rainy 
season (Neckel-Oliveira et al. 2011), when both adult fe¬ 
males and males occupy the lowest strata in the vegetation 
for mating. 
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